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Abstract 

We present an updated version (Nijm93) of the Nijmegen soft-core poten- 
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. tial, which gives a much better description of the np data than the older 

, version (Nijm78). The P^r datum is 1.87. The configuration-space and 

I momentum-space versions of this potential are exactly equivalent; a unique 

0^ ' feature among meson-theoretical potentials. We also present three new AW 

potential models: a non-local Reid-like Nijmegen potential (Nijm I), a local 
version (Nijm II), and an updated regularized version (Reid93) of the Reid 
^ , soft-core potential. These three potentials all have a nearly optimal 

pH I datum and can therefore be considered as alternative partial- wave analyses. 

All potentials contain the proper charge-dependent one-pion-exchange tail. 

^ : 13.75. Cs, 12.40.Qq, 21.30.-hy 
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I. INTRODUCTION 



In the past many nucleon-nucleon (AW) potentials were constructed, which were sup- 
posed to fit the NN scattering data available at the time of construction. The older models, 
from the 1950s and 1960s, are no longer suitable for describing the present set of more nu- 
merous and much more accurate data without refitting the parameters. Out of the various 
potential models constructed in the 1970s, the better ones fitted the data with x^/-^data of 
about 2, where iVdata denotes the number of NN scattering data available at that time in 
the 0-350 MeV energy range. The potentials constructed in the 1980s have only slightly 
improved on this in the sense that, although they have been fitted to try to describe the 
newer and much more accurate data, these models still have x^/^data ~ 2. This number 
should be compared with Xmin/^data = 0.99, obtained in the recently finished Nijmegen AW 
multienergy partial-wave analysis [0] (PWA93) of all pp and np scattering data below 350 
MeV. On statistical grounds, Xmin/^data ~ 1 is about the best one can expect to get in 
partial-wave analyses or for potential models. 

In a recent paper 0, we investigated the quality with respect to the pp scattering data 
below 350 MeV of a number of AW potentials that had appeared in the literature. We found 
that only a few of the potential models we investigated are of a satisfactory quality. These 
models are the Reid soft-core potential Reid68, the Nijmegen soft-core potential |^ 
Nijm78, and the new Bonn pp potential Bonn89. The latter is a readjustment of the 
momentum-space full Bonn potential , in order to fit the pp data. If we do not consider 
the very low-energy (0-2 MeV) pp data, also the parametrized Paris potential Paris80 
gives a satisfactory description of the data. The results of Ref. [@] indicate that, at present, 
the best potential models have x^/^pp ~ 1-9, where Npp denotes the number of pp scattering 
data. Moreover, only models which have explicitly included the pp data in their fit belong to 
this category. Potential models which have been fitted only to the np data often give a poor 
description of the pp data, even after applying the necessary corrections for the Coulomb 
interaction. In Ref. 0] we have demonstrated that most np potentials unfortunately do not 
automatically fit the pp data, a fact which has been generally overlooked. Any WW potential 
should be fitted to the pp data as well as to the np data in order to be able to describe all 
WW scattering data. 

Over the last decade the quality of the np data has increased considerably. Consequently, 
the older potentials (Reid68, Nijm78, Paris80) do not fit these data very well. Also, the much 
newer Bonn potentials already needed revisions and updates [^H. In this paper we present 
updates of the Nijm78 and Reid68 potentials, denoted by Nijm93 and Reid93, respectively. 
Because our analysis of the np data (and hence our careful scrutiny of the np data) has only 
recently been finished we originally constructed an update (Nijm92pp) of the Nijm78 
potential for the pp data only. This pp potential was used in our earlier preliminary np 
analyses P|,|TD|,[TT]] to parametrize the isovector partial waves. It has x^/^pp = 1-4, which 



is not as good as the Nijmegen PWA93. One can wonder whether it is at all possible to 
construct a new class of potential models which fit the WW data with the almost perfect 
X^/Wdata ~ 1- The answer turns out to be affirmative. This could already be surmised from 
the Nijmegen PWA93, because this analysis is in essence an energy-dependent potential 
fitted to the scattering data. (The reason for us using an energy-dependent potential is 
nothing more than just convenience.) In the partial-wave analysis we need 39 parameters to 
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reach x^/^data = 0.99, whereas a conventional potential model typically has only 10-15 free 
parameters. It is therefore perhaps not surprising that the 15-parameter update (Nijm93) 
of the Nijm78 potential, which fits the AW data with x^/^data = 1-87, cannot compete in 
quality with the Nijmegen PWA93. To obtain a high-quality potential we decided some 
years ago to follow a different approach. 

Because the Nijm92pp potential already gives a reasonable description of the pp data, this 
model forms the basis for the construction of a high-quality potential, which can compete 
with the Nijmegen PWA93. We adjust in each partial wave separately only a few of the 
parameters of this potential |]T2|. This way we will be able to construct a potential model 
which fits the data with x^/^data ~ 1- The resulting Reid-like potential Nijm I gives a very 
good fit to the data with x^/^data = 1.03. 

The Nijm I potential contains momentum- dependent terms (as do the Nijm78 and Nijm93 
potentials), which in configuration space give rise to a non-local structure {A(p{r)+(f{r)A) to 
the potential. We also constructed a purely local Nijm II potential, where these momentum- 
dependent terms were intentionally omitted. This local potential Nijm II gives an equally 
good fit to the data as the non-local potential Nijm I. Finally, we constructed a regularized 
update of the Reid68 potential 0, called Reid93. This Reid93 model is also a local potential 
and fits the scattering data very well. These latter three potential models are in a sense 
also alternative partial-wave analyses, because they have roughly the same number of fit 
parameters as our Nijmegen PWA93, these parameters were fitted to the same database, 
and the potential models achieve nearly the same values of Xmm ^^e Nijmegen PWA93 (i.e., 
close to the expectation value). Hence, the differences between, e.g., the phase parameters 
of these models provide an indication for the systematic error in the Nijmegen partial-wave 
analyses. 

In Sec. y we briefly discuss some general features of AW potentials. In Sec. |ITI| we give 
more details regarding the explicit form of the potentials used in this work. Two of these 
potentials are based on the original Nijm78 potential, whereas the third is a regularized 
update of the Reid68 potential in the sense that also in this new model each partial wave is 
parametrized by a number of Yukawa functions. In Sec. IV we discuss the fitting procedure 
and the potentials are presented in more detail. 



II. GENERAL OUTLINE 

The NN potential can be described in momentum space and in configuration space. Since 
it is difficult to solve the full four- dimensional scattering equation, it has become common 
practice first to make a reduction to a three-dimensional scattering equation. Various choices 
are possible, and it is important to note that the potential derived within the chosen reduc- 
tion scheme should only be used in the scattering equation corresponding to that particular 
reduction scheme. These three-dimensional scattering equations can always be written in 
the form of the momentum-space version of the Lippmann-Schwinger equation. If the kine- 
matics is treated relativistically, this is called the relativistic Lippmann-Schwinger equation. 
In configuration space, the differential form of this integral equation is the Schrodinger 
equation. 

The configuration-space potentials are to be used either in the nonrelativistic or the 
relativistic Schrodinger equation 
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{A + k'^)iP = 2MrV ip , (1) 

where A is the Laplacian, and where (non)relativistic refers to the kinematics. For non- 
relativistic kinematics the relation between the center-of-mass energy E and the center-of- 
mass momentum squared k'^ reads E = k'^/2Mr, whereas for relativistic kinematics it reads 



E = ^k^ + Ml + + M| - Ml - M2. 

The earhest potential models were configuration-space potentials to be used in the non- 
relativistic Schrodinger equation. They were phenomenological or semiphenomenological 
parameterizations, based on a general form for the potential. The potential must be invari- 
ant under rotations, reflections, and time reversal, and can be written [|1^] as the sum of 6 
independent terms, V = ViPi. A common choice for the 6 operators Pi in configuration 
space is 

P2 = (71^2, 

P3 = S12 = S{ai- r) (0-2 ■ r) - (a^ ■ 0-2) , 
P4 = L ■ S, 

= Q12 = i[(o-rL)(o-2-L) + (cr2-L)(<Ti-L)], 

P6 = i(o-i-<T2)-L. 

These operators are also frequently referred to as the central, spin-spin, tensor, spin-orbit, 
quadratic spin-orbit, and antisymmetric spin-orbit operators, respectively. For identical- 
particle scattering, the antisymmetric spin-orbit operator Pg cannot contribute, whereas Vg 
vanishes when charge independence is assumed (which is usually the case for AW potential 
models). In general [0, each potential form Vi in configuration space is a function of 



and of the operators and L^. In most approaches one only keeps the dependence on r^, 
while the dependence (when included) is often only present in a linear way in the central 
potential Vi. The inclusion of the Q12 operator was found to be necessary, because otherwise 
it was impossible to describe simultaneously the ^5*0 and ^D2 phase shifts using the same 
static potential. The presence of the operator Q12 in the potential can to a certain extent 
be simulated by introducing non-local potentials ||T4 |. 



In the expansion I]i=i ViPi, the potential forms Vi are generally assumed to be the same 
in all partial waves. The potential differences between the partial waves are dictated by the 
differences in the expectation values of the operators Pj in these partial waves. The Reid68 
potential , however, is based on a quite different approach. Rather than having 6 potential 
forms Vi which are the same for all partial waves, now each partial wave is parametrized 
separately. The potential forms Vi therefore not only depend on and L^, but also on S'^ 
and J^. In this paper we present some potentials based on this approach that each partial 
wave is parametrized independently. We refer to these models as Reid-like models. 

With the discovery of the heavy mesons in the 1960s, it became common practice to 
write the potential as a sum over one-boson-exchange (OBE) potentials. The expressions 
for these OBE potentials are usually derived in momentum space. Introducing momentum 
vectors 

k = P/-Pi, q=|(P/ + Pi), n = qxk. 
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in terms of initial (pi) and final (p/) momenta, the equivalent in momentum space to Eq. 
reads 



p _ 
1 — 


1 
1 
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|(cri - (T2 


)-n. 



(3) 



The potential forms Vi in momentum space are functions of k, q, n, and the energy. Although 
Eq. (0) provides an adequate set of 6 linearly independent operators, the Qu operator in 
configuration space is not the exact Fourier transform of the {cti ■ n) (cr2 ■ n) operator in 
momentum space. This is of importance if we want both the momentum-space and the 
configuration-space versions to produce exactly the same phase shifts and bound states, 
which is only possible when the configuration-space version is the exact Fourier transform of 
the momentum-space version, and vice versa. This implies that we have to use the inverse 



Fourier transform of the Q12 operator; i.e., the potential contribution (cti ■ n) (cr2 ■ n) V5(k^) 
is to be replaced by 

P,V,ik') - pU"^' dk''V,{k'') , (4) 

J 00 

where 

P5 = [(crrq)(cr2-q) - q^(CTrcr2)] 

-i[(^l-k)(^2-k)-k2(^l-a2)] . (5) 

Other restrictions imposed on the momentum-space potential forms Vi in that case are that 
they should not depend on the energy, while the q dependence should be of second order at 
most (see also below). 

When the potentials are evaluated in momentum space and then Fourier transformed 
to configuration space, they are usually first regularized to remove the singularities at the 
origin. This can be achieved by introducing a form factor F(k^). A typical Fourier transform, 
encountered in transforming the momentum-space potential to configuration space, then 
reads 



. (6) 



m 
An 



The results for various frequently used choices are: 

(i) No form factor at all, -F(k^) = 1- This yields the familiar Yukawa function 



r] 



e-'"7mr , (7) 
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and the singularities at the origin are still present; 

(ii) Monopole form factor, -F(k^) = (A^ — iv?)/{h? + k^), normalized such that at the pole 



-m 



1. This yields 



/mr ; 



(iii) Dipole form factor, -F(k^) = (A^ — m^)^/(A^ + k^)^, yielding 



A^ 



m 



(iv) Exponential form factor, F(k^ 



e-'""erfc 



-e™"erfc 



2A2 

yielding 

m Ar\ 

A ~ tJ 

m kr 



-Ar 



mr 



(8) 



(9) 



)1/ 



2mr 



(10) 



where erfc(x) is the complementary error function 

2 



erfc(x) 



dt e 



We follow the normalization of Ref. Q]. This means that for the exponential form factor 
F(0) = 1. 

Because k^ can be written as (k^ + m^) — m^, we find that in the absence of a form factor 



0^(r) = 0^(r) - 47r(5^(mr) 



When there is a form factor, this relation still holds, but the ^-function contribution is 
smeared out. 

Using our definition (|^), the Fourier transforms for the tensor and spin-orbit potentials 
can be simply expressed in terms of derivatives of the central function, i.e.. 



0O(r) 



3m^ dr ( r c/r ) ^'"^ ^ 



(12) 



r dr 

In order to ensure regularity at the origin for the tensor and spin-orbit functions, one must 
choose at least the dipole or exponential form factor. In that case, the tensor function also 
vanishes at the origin, as it should. 

The presence of explicit momentum-dependent terms in the momentum-space potential 
gives rise to non-local structures in the potential in configuration space. The terms pose 
no difficulties for the configuration-space potential as long as they are linear in q^. The 
typical Fourier transform of such a term is given by 



d^k 



ikr 



(27r)3 k2 + m2 



+ ik^ i^(k^ 



m 
47r 



2 Mr 2Mr 



(13) 
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where (/9(r) = M,.0^(r). It is well known how to handle such a (Ayj + yjA) term [[16 
The absence of terms in the momentum-space potential will result in a radially local 
configuration-space potential. 

The three new potential models (Nijm93, Nijm I, and Nijm II) presented in this paper 
are based on the original Nijm78 potential with the exponential form factor, whereas the 
update of the Reid68 potential (Reid93) is regularized using a dipole form factor. 



III. STRUCTURE OF THE POTENTIALS 
A. One-pion-exchange potential 

An important feature of the potential models presented in this paper is that in the 
one-pion-exchange (OPE) part of the potential, we explicitly distinguish between neutral- 
pion and charged-pion exchange. The pion masses are m^o = 134.9739 MeV and 
mjr± = 139.5675 MeV. Almost all other potentials that have appeared in the literature use 
a mean pion mass. In these other models the isovector np phase parameters are larger in 
magnitude than the corresponding pp phase parameters. By explicitly including the pion- 
mass differences exactly the opposite occurs: the isovector np phase parameters are smaller 
than the corresponding pp phase parameters. This is a unique feature of the potentials 
presented here. 

Defining 

V{m) = i m[0°(m,r)^i2 + |0j;(m,r)(<Trcr2)J , (14) 

the OPE potential for pp scattering is given by 

Vope(pp) = /.V(m.o) , (15) 

whereas for np scattering it reads 

VoMnp) = -/.V(m^o) ± 2f^V{m^±) , (16) 

where the plus (minus) sign corresponds to total isospin / = 1 (0). The scaling mass m^ri in 
V{m) is introduced in order to make the pseudovector coupling constant dimensionless. 
It is conventionally chosen to be equal to the charged-pion mass. The explicit distinction 
between neutral and charged pions implies that the isovector np and pp OPE potentials are 
different, and so charge independence is broken. In our present models we assume, however, 
that the pion-nucleon coupling constants obey charge independence. 



B. Nijmegen potential 

In this section we briefly discuss the structure of the Nijmegen potential. More details can 
be found in Refs. [§,0. The basic functions are the one-boson-exchange (OBE) potential 
functions with momentum-dependent central terms and exponential form factors. The meson 
exchanges we include are those due to pseudoscalar mesons (tt, rj, rj'), vector mesons (p, u, (p), 
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and scalar mesons (oq, /o, e). Here we use the modern nomenclature for the scalar mesons, 
i.e., ao(983) corresponds to the 6 of Ref. [Q, and /o(975) to the S*. The e meson would 
correspond to an /o(760). No such meson is listed by the Particle Data Group |T^; however. 



a recent analysis of the vrA^ — * %^n~N reaction [^] provides evidence for a scalar-isoscalar 
resonant state 0'*'^(750). In the Nijmegen potentials the e meson corresponds to a broad 
meson (see below) where the pole in its propagator is chosen to correspond to the pole 



position in the complex energy plane of the isoscalar tttt 5* wave |^ . Here we will retain the 
name of e meson. The aforementioned meson exchanges can be identified with the dominant 
parts of the lowest-lying meson trajectories in the complex J plane. We furthermore include 
the dominant J = parts of the Pomeron, and of the /2, f'21 and 02 tensor-meson trajectories. 
They give rise to Gaussian potentials. 

The meson propagators including the exponential form factor read 

A(k^m^A^) = -^e-'^VA^ (17) 
For the Pomeron-type exchanges we have 

A{k^ml) = ^e-^y'-'^, (18) 
p 

where rup has the dimension of a mass and will be called the Pomeron mass, and Mp is 
a scaling mass, chosen to be the proton mass. The different potential forms are evaluated 
in momentum space and the resulting expressions are essentially those of Refs. 0] (save 
some misprints [^) with the following differences: (i) We explicitly account for the proton 
and neutron mass difference; (ii) the differences between the neutral and charged pion (see 
Sec. [HI A|) , and between the neutral and charged p meson are explicitly included; (iii) we have 
adjusted the quadratic spin-orbit operator of the potential in momentum space to include 
the P5 contribution as in Eq. (|^). The effect of the first modification is obviously rather 
small. The second modification (as well as the first) implies that charge independence is 
broken in the non-OPE part of the potential as well. For pp scattering the potential consists 
of only neutral-meson exchange, Vpp = V^(neutral), whereas for np scattering it consists of 
neutral- meson and charged- meson exchange, depending on the total isospin as in Eq. (p!6[), 
so Vnp = — \^(neutral) ± 2l^(charged). This distinction replaces the factor (ti-T2) used in 
the old Nijm78 potential. Finally, the third modification means that we have constructed 
a potential which is exactly equivalent in both momentum space and configuration space, 
a unique feature of these Nijmegen potentials. For example, for the parametrized Paris 
potential [0 this is not true since it uses the same parameters in combination with the Q12 
operator in configuration space as it does with the {cri ■ n) {(T2 ■ n) operator in momentum 
space. 

Next we briefly discuss the coupling constants. For definitions and references we again 
refer to Refs. P,pr8[. The coupling constants of the pseudoscalar mesons are related via 
SU(3) and singlet-octet mixing. The octet coupling /^^ is calculated using ap = 0.355 
for the F/{F + D) ratio. For the singlet-octet mixing angle we use dp = —23° to define 
the physical coupling constants frj and /^/. This leaves the singlet coupling /^^ and the 
pion coupling as free input parameters. However, in our partial-wave analysis of the pp 
scattering data [^, we found for the ppn^ coupling constant = 0.0749(7). This value 
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was later confirmed in a combined partial-wave analysis of all pp and np scattering data, 
assuming charge independence for the pion-nucleon coupling constants There the 

value = 0.075 is recommended for the pion-nucleon coupling constant at the pion pole. 
This is the value we adopt in our construction of the new Nijmegen potentials, and so it is 
not included clS 8b free parameter. 

For the vector mesons we assume that the p meson is universally coupled to the isospin 
current {ay = 1) to define the octet coupling gvg. For the singlet-octet mixing angle we take 
9v = 37.5°, which fixes the physical coupling constants and g^p in terms of gp and gy^- 
The (p meson is assumed to have = 0. The free parameters are now gp, gv^, fp, and 

For the scalar mesons we do not apply any constraints for the coupling constants, since 
the singlet-octet mixing angle for the scalar mesons is still an unsettled problem (see also 
Ref. [|18|). The free parameters are the oq, /o, and e coupling constants. 

For simplicity we take a single mass parameter rrip for the Pomeron, and for the J = 
parts of the /2, /a, and 02 tensor-meson trajectories. We use two coupling constants: ga2 for 
the isovector 02 meson and gp for the isoscalar and Pomeron exchanges. 

For each type of exchange we use an independent cutoff mass, so we have three cutoff 
parameters Ap, Ay, and A5. This brings us to a total of 14 free parameters. 

We conclude this section with a discussion of the treatment of the broad p and e mesons. 
The width of a broad meson can be accounted for [^,^,^ by replacing the propagator 
A(k^) = l/(k^ + m^) of a stable meson by a dispersion integral 

A(k=)=r4?^. (19) 



with mass distribution 

P^"^^ = — T^lfi ' (20) 



7 = mr(m2-m2)-("+i/2) . (21) 



and where 



Here F denotes the width and n = 0, 1 for spin-0 and spin-1 mesons, respectively p3[. 
The charged p meson decays into a neutral and a charged pion and the threshold mass is 
rrit = m^o + m^ri . The neutral p meson cannot decay into two neutral pions and it decays 
into two charged pions, and so now the threshold mass is = 2m^±. The e meson is an 
isoscalar meson which decays into both two neutral or two charged pions in the ratio 1:2. In 
our present models these distinctions have been explicitly accounted for, which is another 
extension of the old Nijm78 model. 

The configuration-space potential due to the exchange of a broad meson is calculated 
exactly. This exact potential is then approximated by the sum of two potentials of stable 
mesons [^1 



mt 



dm'2m' p[m''^)m' (lP(j{m' , r) 

<P^c{mur) + (32m2<PU^2,r) . (22) 
Fitting from 0-2 fm yields the values as given in Table |. 
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C. Regularized Reid potential 



A disadvantage of the original Reid68 potential is that, at the time of its construction, 
the quality of the np data was very poor. As a consequence, the Reid68 potential can no 
longer properly describe the numerous new and much more accurate np data. Another 
disadvantage is that the Reid68 potential has an singularity in all partial waves. Here 
we present an updated version of the Reid potential, where these singularities have been 
removed via the inclusion of a dipole form factor. With this choice, the tensor potential now 
also vanishes at the origin, as it should. 

As is the case for the original ReidGS potential, the OPE potential is explicitly included, 
while we now account for the neutral-pion and charged-pion mass differences as in Eqs. (p!^- 
(p!6|). For the pion-nucleon coupling constant at the pion pole we take P,[TT| fl = 0.075, 



and for the dipole cutoff parameter we choose A = 8m^. In the OPE potential (JT^) we use 
0^ only for the 5* waves. For all other partial waves, we found it more convenient to use 
0^ instead of (f)}j. Note that 0^ equals ipQ up to a modified 6 function [see Eq. (pH])], and 
that this modified 6 function is screened by the centrifugal barrier for all these other partial 
waves, except the S waves. 

Starting with this OPE potential, the potential in each partial wave can now be extended 
by choosing a convenient combination of central, tensor, and spin-orbit functions with arbi- 
trary masses and cutoff parameters. In the construction presented in the following we (more 
or less arbitrarily) settled for integer multiples of a mean pion mass m = (m^o + 2m^±)/3, 
while the cutoff mass in the dipole form factor is chosen to be A = 8m everywhere. For 
notational reasons we next define 

Y{p) = prn(l)Q{pfn,r), 
Z{p) = pm0° (pm, r), 
W{p) = prn(f)gQ{pm, r), 

with p an integer and (p^x given by Eqs. and ([l^). For the coefficients multiplying these 
functions, we use Aip for the isovector potentials, whereas the coefficients Bip are for the 
isoscalar and np ^Sq potentials. The index i subsequently labels the different partial waves. 
For the total potential in a particular partial wave one should, of course, add the appropriate 
OPE potential as given by Eqs. (|T^-(|TB]). 



For the non-OPE parts in the isovector singlet partial waves {I = 1, S = 0, L = J) we 



use 



VppCSo) = ^12^(2) + Ai3F(3) + ^14^(4) 

+A,,Y{5) + AieY{Q) , 
VnXSo) = 5i3F(3) + BuYiA) + ^15^(5) + B,,Yi6) , 

VCD,) = ^24^(4) + ^25^(5) + ^26^(6) , (23) 
VCC) = ^33^(3) , 

VCJj) = VppCSo) forJ>6, 

where the distinction between the pp and np ^Sq potentials is necessary because of the well- 
known breaking of charge independence in the pp and np ^Sq partial waves. The coefficients 
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Aip and Bip are to be fitted. The presence of tlie two-pion range piece Aq2Y{2) in tlie pp 
^So potential is purely coincidental, and was only included to improve the quality of the fit. 
A similar term in the np ^Sq was much less effective, and so we decided to leave it out. 
For the non-OPE parts in the isoscalar singlet partial waves {I = 0, S = 0, L = J) we 

use 

VCPi) = B23Y{3) + 524^(4) + 525^(5) + 526^(6) , 

VCF,) = Bs3Y{3) + Bs,Y{5) , (24) 
VCJj) = VCPi) for J > 5 . 

For the isovector triplet uncoupled partial waves {I = 1, S = 1, L = J) we use 

VCPo) = A,sY{3) + ^45^(5) + ^4.3^(3) , 

VCP,) = A,sY{3) + A,,Y{5) + A,,sZ{3) , (25) 
VCF,) = A,,Yi3) , 

and the isoscalar triplet uncoupled partial waves {I = 0, S = 1, L = J) are parametrized as 

VCD2) = B,sY{3) + B,,Y{5) + 54,3^(3) , 

VCC) = ^53^(3) . (26) 

Following the parametrization of the original Reid68 potential, the non-OPE potential in 
the triplet coupled partial waves (S" = 1, L = J ± 1) is parametrized as 

V = Vc + VtS^2 + V^soL ■ S , (27) 

where the isovector (/ = 1) potentials are given by 

Vc = Ar,Y{3) + Ar,Y{A) + Aj,Y{5) + Ar,Y{6) , 

Vt = Ar,,Z{A) + AueZiQ) , 
Vso = Ar^sWi3) + Aj^,Wi5) for J = 2 , (28) 
Vso = As^sWi3) forJ = 4, 

and the isoscalar (/ = 0) potentials read 

Vc = Be2Y{2) + 563^(3) + ^64^(4) 

+^65^(5) + ^66^(6) , 

Vt = Be,,Z{A) + B^,qZ{Q) , (29) 
Vso = B^^^W{3) + 56^5W^(5) for J = 1 , 
Vso = B,^^W{3) + B,^,W{^) for J = 3 . 

Finally, for the triplet isovector partial waves (/ = 1, S* = 1) with J > 5 we use Eq. ( pTj) 
with the central and tensor potentials of Eq. (pSf), and the spin-orbit potential equal to zero. 
Similarly, for the triplet isoscalar partial waves (/ = 0, 5* = 1) with J > 5 we use the central 
and tensor potentials of Eq. (|29|). This choice is analogous to the extension of the Reid68 
potential to the higher partial waves as given by Day . 
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IV. RESULTS 



The parameters of the potential models are optimized by minimizing the in a direct 
fit to the data. Since the scattering data are spread over a large number of energies (about 
200 different energies for the data and almost 400 different energies for the np data in 
the 0-350 MeV energy range) and the phase parameters need to be calculated up to at 
least J ~ 6, the Schrodinger equation then has to be solved a very large number of times. 
This approach, therefore, is not very practical when the model parameters do not yet have 
reasonable values and, consequently, the is still very high. A more convenient approach is 
to start with the Nijmegen representation |I| of the hypersurface of the scattering data. 
It is obtained from the 10 single-energy analyses and consists of 10 sets of phase parameters 
and the error matrix, each at a different energy. The error matrix is the inverse of half the 
second-derivative matrix of the hypersurface with respect to the phase parameters up 
to J = 4 within the energy bin of the single-energy analysis. This hypersurface is, in 
principle, independent of the particular partial-wave analysis. In practice the representation 
we use is somewhat dependent on the Nijmegen multienergy analysis. The crucial point 
is, however, that it provides a very good and concise representation of the scattering data. 
For each change in the model parameters we need to solve the Schrodinger equation for all 
partial waves up to J = 4 at only 10 different energies, which allows for a much quicker 
optimization for the parameters of the potential model. 

In the last stage of the fitting procedure the potential parameters have been further 
optimized in a very time-consuming direct fit to the data. In this case we use the potential 
model in all partial waves with J < 6, whereas in the higher partial waves we include OPE 
only. The final Xmin of ^^e potential should be obtained from this direct comparison with 
the experimental data. 

A. Nijm92pp 

Our first improvement of the Nijm78 potential was already started several years ago, 
when we constructed an update to the data of the Nijm78 potential. This potential 
has been used in the Nijmegen analyses |^Jl^|Tl| to parametrize the np isovector lower 
partial waves (except the np ^S'o). In our latest analysis [Q (PWA93) we refer to it as the 
Nijm92pp potential. We found that a good fit to the pp data could be obtained using one 
cutoff parameter A = 827.53 MeV for all three types of meson exchanges. Some of the 
coupling constants were not refitted, but were kept at the values of the original Nijm78 
potential. The reason is that, when we only fit to the pp data, we cannot incorporate the 
isospin dependence of the isovector exchanges (there are only 7 = 1 partial waves). A direct 
comparison of this Nijm92pp potential with the pp data yields x^ = 2487.1 for 1787 data, 
which means I^pp = 1-4- 

B. Nijm93 

Now that the Nijmegen analysis of the np data is also finished ^ and we have carefully 
scrutinized the np database, we can extend the update of the Nijm78 potential to include 
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the fit to the np scattering data as well. This model we refer to as Nijm93. 

As already mentioned in Sec. Ill C , the np ^Sq partial wave has to be parametrized 
separately. The reason is that there is clear evidence for breaking of charge independence 
in the ^5*0 scattering lengths ttpp and a„p. This difference in scattering lengths carries over 
into an approximately 2° phase-shift difference between the pp and np ^Sq phase shifts at 
higher energies. This difference cannot be explained as being only due to the difference 
between the pp and np OPE potentials. Allowing for a different value for the neutral-pion 
and charged-pion coupling constants does not help either, because the scattering lengths are 
very insensitive to these kind of changes. To accommodate the pp and np ^Sq differences, 
we therefore introduce a purely phenomenological breaking of charge independence between 
the p° and coupling constants. This breaking of charge independence is only assumed in 
the ^So partial wave; for all other partial waves the p° and coupling constants are taken 
to be the same. 

The parameters for the Nijm93 potential, rounded to four or five significant figures, can 
be found in Table 0, where the meson masses are the masses as listed by the 1990 Particle 
Data Group [|T^. The coupling constants are the values at = 0. The pion coupling 
constants are fixed at = 0.075 at the corresponding pion pole = — m^o or k^ = — m^±; 
hence the different entries for neutral- and charged-pion exchange at k^ = 0. The e and 
Pomeron coupling constants are rather large, whereas the u coupling constant is reasonably 
small. For the p coupling constants we find {f/g)p = 4.094, which is close to the value 
3.7 from naive vector-meson dominance of the isovector electromagnetic form factors of the 
nucleon. 

For the Nijm93 potential we find x^ipp) = 3175.6 for 1787 pp data and x^(np) = 4848.4 
for 2514 np data. So for the pp data x^/^pp = 1-8) ^oi the np data x^/^np = 1-9, and 
for all AW data x^/^data = 1-87. We find that this 15-parameter conventional meson- 
exchange potential cannot do better than x^/^data ~ 1.87, a result which is also found 
for similar potential models such as the Paris80 and all the Bonn potentials. Apparently, 
the conventional meson-exchange potentials cannot compete in quality with the Nijmegen 
PWA93. This indicates that these models lack some important physics. 



C. Nijm I and Nijm II 

In order to be able to construct a potential model which is of almost the same quality 
as the Nijmegen PWA93 (x^/^data ~ 1), we follow a different approach and take advantage 
of the success of the Reid68 potential. We expect that, when we start with the Nijmegen 
potential (which already has a reasonable x^ on the pp data), we can construct a Reid-like 
potential where in each partial wave we probably need to adjust only a few parameters in 
order to arrive at x^/^data ~ 1- The potential forms are then given by a set of slightly 
adjusted Nijmegen potentials, each representing one particular partial wave. Starting with 
the parameters of the Nijm92pp potential [jl2|, we find that for most partial waves an ad- 



justment of the fp and coupling constants already gives very good results. Counting all 
parameters which have been adjusted in the fit of each partial wave, we arrive at a total 
of 41 parameters. This should be compared with the 39 parameters used in the Nijmegen 
PWA93. 
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In the last stage, the parameters of this Reid-hke potential are optimized in a direct fit 
to the data. The potential is then used in all partial waves up to J = 6 simultaneously. 
This model we refer to as Nijm I. It has x^iPP) = 1795.8 and x^(^p) = 2627.3, and so 
X^/iVdata = 1-03 on all pp and np scattering data. 

We have also constructed a local Reid-like Nijmegen potential, where we leave out 
the explicit momentum-dependent terms which give rise to non-local contributions to the 
configuration-space potential as expressed in Eq. ([TB|). We follow the same procedure as 
for the non-local Nijm I potential. First, the parameters are adjusted in a fit to the rep- 
resentation of the hypersurface, and then further optimized in a direct fit to the data. 
For this local potential, denoted by Nijm II, we use a total of 47 parameters and we find 
X^ipp) = 1795.8 and x\np) = 2625.7, and so xVA^'data = 1-03. 

Although these potentials are purely phenomenological (except for the correct OPE tail) 
and the coupling constants have no physical meaning, these potentials are the first to give an 
excellent description of the AW scattering data. They have already been used successfully 
in three- and many-body calculations 0)11 • 



D. Reid93 



Finally, we have constructed an updated Reid potential based on the original Reid68 
potential. This regularized version, denoted by Reid93 and discussed in Sec. [Ill C| , gives 



an equally good description of the data as do the Nijm I and Nijm II potentials. The 50 
phenomenological potential parameters Aij and Bij were fitted to the data, resulting in 
X^{pp) = 1795.1 and x^(^p) = 2624.6, and so also for this potential x^/^data = 1-03. This 
Reid93 potential has been used in a triton calculation as well 



E. Comparison of the potentials 



The results of the non-local Nijm I potential, the local Nijm II potential, and the local 
Reid93 potential are summarized in Table |T|. The results of the Nijmegen PWA93 are 
shown for comparison. Although x^/^data for these three potentials is already very good, 
their description of the np data is still not as good as the description of these np data 
in the partial-wave analysis. Here we should also mention that we did not do a thorough 
investigation into the minimum number of parameters required to get these results (as we did 
for the partial- wave analysis). The reason is that in order to do this properly, all potential 
parameters have to be fitted simultaneously to all NN data. But we found it more successful 
to do a large number of fits where in each separate run only a (completely arbitrary) subset 
of the potential parameters was optimized. Therefore, we cannot rule out the possibility 
that an equally good fit can be obtained with a few parameters less. 

As already mentioned in the Introduction, these new potentials (except Nijm93) are in 
a sense alternative partial-wave analyses. The differences between the phase parameters of 
the potentials and the phase parameters of the Nijmegen PWA93 are shown in Tables ^ 
and 0. These differences provide an indication for the systematic error on the results of the 
Nijmegen PWA93. For the np phase parameters the differences at some energies (especially 
with the Reid93 potential) are relatively large. However, one has to bear in mind that x^{f^p) 
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of the potential models is still substantially higher than that of the multienergy partial- wave 
analysis. On the other hand, the variation in the mixing parameter ei is small. It has often 
been claimed that this mixing parameter is ill determined and a wide range of values from 
potential models seemed acceptable (see, e.g., Ref. [^H). However, in our publication of 
the Nijmegen PWA93 we already argued that ei is in fact known very accurately. This is 
confirmed in Fig. |I|, where we note that the results of the Nijm I and Nijm II models lie 
essentially within the statistical uncertainty as obtained in the Nijmegen PWA93. Above 
150 MeV, the result of the Reid93 model rises too strongly but is still within 2.5 standard 
deviations of the Nijmegen PWA93. 

The potentials between and 2 fm for the singlet, triplet uncoupled, and triplet coupled 
np partial waves are shown in Figs. ^, ^, and ^ respectively. For the non-local potential 
Nijm I we plot V/{1 + 2(f), which more or less represents the effective potential when non- 
local terms are present (see Refs. |T^J^). For coupled channels, the potential is a 2 x 2 



matrix and the ei and €2 plots in Fig. ^ represent the off-diagonal elements of the potential. 
The main differences between the potentials show up in the inner region, i.e., for r < 1 fm. 
In general, the non-local Nijm I potential is much softer than the local Nijm II and Reid93 
potentials, while the Nijm II potential is again much softer than the Reid93 potential. The 
reason for this softness of the Nijmegen potentials is the exponential form factor. 

Finally, all potential models have been fitted to the deuteron binding energy B = 
2.224 575(9) MeV p9| using relativistic kinematics, i.e.. 



rather than B = /2Mr. We also constructed versions of the Nijm I and Nijm II potentials 
to accommodate the latter nonrelativistic form. In any case the value B = 2.224 575 MeV 
is exactly reproduced to this accuracy. Some of the other deuteron parameters are listed in 
Table |V|. The different potential models all give very similar results. Because we consider 
the potentials Nijm I, Nijm II, and Reid93 as alternative partial wave analyses, the values 
of the deuteron parameter 77, As, and A^^ given by these potentials are new experimental 
determinations of these quantities. For the D/S-state ratio 77 we find 77 = 0.0252(1) in good 
agreement with the recent determination by Rodning and Knutson |^0|] of 77 = 0.0256(4). 
For the asymptotic S'-state normalization A 5 we obtain As = 0.8843(10) fm~^/^, which is in 
agreement with the determination by Kermode et al. of As = 0.8883(44) fm~^/^. This 



results in A^^ = AKl + r]"^) = 0.7825(20) fm~^. However, there have been other experimental 
determinations of these quantities which are not always in agreement with the values quoted 
above. For a more complete list of these experimental determinations and a discussion of 
the differences between them, we refer to Ref. and references cited therein. A direct 
comparison of our value Qd = 0.271(1) fm^ of the deuteron quadrupole moment with the 
experimental value 0.2859(3) fm^ of |]53| is only possible after all possible corrections have 
been accounted for, which is outside the scope of this paper. However, we would like to 
turn the argument around and suggest that these corrections must obviously be about 
0.015 fm^. It is quite interesting to see that for our best potentials the Z)-state probability is 
Pd = 5.665(30)%. The deuteron parameters as well as the results for the scattering lengths 
for both the potential models and the Nijmegen PWA93 will be discussed in more detail 
elsewhere ||3^ . These potentials were used in calculations of the triton binding energy. 
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It turned out that all these potentials underbind the triton by roughly 800 keV, a result 
which can be expected from their values. 



V. CONCLUSIONS 

In this paper we presented an update Nijm93 of the old Nijm78 NN potential. It contains 
the correct OPE tail and has x^/^data = 1-87. Although it cannot compete in quality with 
the Nijmegen partial-wave analysis (a feature apparently all conventional meson-exchange 
potentials suffer from), the description of the np data of the new Nijm93 model is substan- 
tially better than that of the original Nijm78 potential, which was fitted to the old 1969 



Livermore database |35|. Here we would also like to point out that in the Nijm93 model 
we do not include two-meson-exchange contributions such as vrvr and np exchange, and we 
still get a reasonable description of the ^Pi and phase shifts. This in contrast to claims 
made in the literature |]^,|^ that this is impossible. However, this result is obtained at the 
cost of having a rather large value for the pseudoscalar (pion) cutoff mass of Ap = 1177.11 
MeV. 

We have also presented three new high-quality AW potentials. The Nijm I potential 
is a non-local Reid-like potential where each of the lower partial waves up to J = 4 is 
parametrized separately. For the higher partial waves we use the parameters of the Nijm92pp 
potential, which is an update to the pp data of the original Nijm78 potential. The Nijm II po- 
tential is a local Reid-like potential, and does not contain any explicit momentum-dependent 
terms. Both potentials fit the NN scattering data with a nearly optimal Xmin/^data = 1.03. 
A regularized update of the Reid potential, denoted by Reid93, gives the same excellent 
xLJNd.t. = 1.03. 

Computer codes for these Nijmegen potentials Nijm I, Nijm II, and Nijm93, and for the 
regularized Reid93 potential, in configuration space as well as in momentum space, can be 
readily obtained via anonymous FTP from thef-nym.sci.kun.nl. 
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TABLES 



TABLE L Values for the parameters of Eq. (^) of the two-pole approximation for the broad 
e meson and the broad neutral and charged p mesons. Masses m and widths F are in MeV. 



n 1 1 

m 760.0 768.7 768.3 

r 640.0 152.4 149.1 

/3i 0.16900 0.26552 0.38755 

mi 487.818 645.377 674.152 

/?2 0.61302 0.56075 0.45083 

m2 1021.14 878.367 929.974 



TABLE II. Masses and meson-nucleon coupling constants at = for the Nijm93 potential. 
For the np ^Sq partial wave the coupling constants of the charged-rho meson are increased by 
4.371% (see text). P in the last line denotes the Pomeron. Note that {f /g)p = 4.094. 

m (MeV) f 

139.5675 0.07395 

TfO 134.9739 0.07402 

r] 548.8 0.01514 

r]' 957.5 0.01466 
Ap 1177.11 

p± 768.3, F = 149.1 0.8481 14.217 

/9° 768.7, F = 152.4 0.8481 14.217 

uj 781.95 9.1765 0.3383 

(f) 1019.412 0.0985 
Av 904.50 

a^,aS] 983.3 1.9174 

e 760.0, F = 640.0 28.196 

/o 975.6 12.142 
As 554.40 

02 208.16 0.0486 

P,/2,/2 208.16 27.339 
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TABLE III. foi' the new potential models in comparison with the Nijmegen multienergy 
analysis PWA93. We also show the number of parameters (A'par) for each model. 



PWA93 Nijm I 

pp 1787.0 1795.8 

np 2489.2 2627.3 

Total 4276.2 4423.1 

iVpar 39 41 

xV^^data 0.99 1.03 



Nijm II Reid93 Nijm93 

1795.8 1795.1 3175.6 

2625.7 2624.6 4848.4 

4421.5 4419.7 8023.9 

47 50 15 

1.03 1.03 1.87 



TABLE IV. pp phase 


shifts in 


degrees. For 


each energy the 


rows give the values from the 


AT' • 
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multienergy partial-wave 


analysis Q, 


the non-local Nijm I potential. 


the local Nij 


TT 

m 11 


potential, 


and the Reid93 potential, respectively. 
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0.10 




48.68 


0.70 


8.60 


—4.91 


2.50 


-0.81 


0.11 




48.72 


0.70 


8.52 


—4.89 


2.49 


-0.81 


0.11 




48.71 


0.69 


8.60 


— 4.yu 


z.4y 


-0.80 


0.10 


50 


38.92 


1.71 


11.47 


—8.32 


5.85 


-1.71 


0.34 




38.91 


1.71 


11.55 


-8.31 


5.85 


-1.70 


0.34 




38.92 


1.70 


11.48 


-8.30 


5.84 


-1.70 


0.34 




39.03 


1.68 


11.67 


-8.30 


5.83 


-1.69 


0.34 


100 


24.98 


3.79 


9.45 


-13.26 


11.01 


-2.66 


0.82 




24.96 


3.73 


9.50 


-13.30 


10.96 


-2.63 


0.82 




24.91 


3.75 


9.55 


-13.33 


10.97 


-2.64 


0.83 




25.09 


3.71 


9.79 


-13.30 


10.97 


-2.61 


0.81 


150 


14.77 


5.61 


4.74 


-17.43 


13.98 


-2.87 


1.20 




14.79 


5.60 


4.63 


-17.51 


13.94 


-2.87 


1.19 




14.70 


5.61 


4.77 


-17.54 


13.95 


-2.87 


1.20 




14.83 


5.55 


4.97 


-17.49 


13.96 


-2.85 


1.16 


200 


6.57 


7.06 


-0.37 


-21.25 


15.63 


-2.76 


1.42 




6.66 


7.20 


-0.63 


-21.32 


15.65 


-2.80 


1.39 




6.56 


7.15 


-0.47 


-21.29 


15.63 


-2.79 


1.42 




6.62 


7.08 


-0.32 


-21.26 


15.63 


-2.77 


1.36 


250 


-0.30 


8.27 


-5.43 


-24.77 


16.59 


-2.54 


1.47 




-0.20 


8.50 


-5.72 


-24.81 


16.65 


-2.62 


1.39 




-0.25 


8.41 


-5.61 


-24.67 


16.61 


-2.61 


1.44 




-0.23 


8.35 
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16.59 


-2.55 


1.39 


300 
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17.17 


-2.34 


1.34 




-6.18 


9.52 


-10.49 


-28.02 


17.26 


-2.41 


1.20 




-6.12 


9.42 


-10.49 


-27.71 


17.21 


-2.41 


1.29 




-6.10 


9.40 


-10.29 


-27.86 


17.15 


-2.26 


1.25 


350 


-11.11 


10.69 


-15.30 


-30.89 


17.54 


-2.21 


1.04 




-11.51 


10.28 


-14.94 


-30.98 


17.63 


-2.23 


0.86 




-11.28 


10.24 


-15.08 


-30.45 


17.61 


-2.23 


0.98 




-11.22 


10.30 


-14.80 


-30.77 


17.49 


-1.95 


0.95 
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TABLE V. np phase shifts in degrees. For each energy the rows give the values from the 
Nijmegen multienergy partial- wave analysis the non-local Nijm I potential, the local Nijm II 
potential, and the Reid93 potential, respectively. 











3Pi 


35i 


ei 




^D2 




«P2 


£2 


^P2 


1 


62.07 


0.18 


-0.19 


-0.11 


147.75 


0.11 


-0.01 


0.00 


0.01 


0.02 


-0.00 


0.00 




62.11 


0.18 


-0.19 


-0.11 


147.76 


0.10 


-0.01 


0.00 


0.01 


0.02 


-0.00 


0.00 




62.09 


0.18 


-0.19 


-0.11 


147.75 


0.10 


-0.01 


0.00 


0.01 


0.02 


-0.00 


0.00 




61.89 


0.18 


-0.19 


-0.11 


147.73 


0.10 


-0.01 


0.00 


0.01 


0.02 


-0.00 


0.00 


5 


63.63 


1.63 


-1.49 


-0.94 


118.18 


0.67 


-0.18 


0.04 


0.22 


0.25 


-0.05 


0.00 




63.74 


1.62 


-1.50 


-0.93 


118.19 


0.67 


-0.18 


0.04 


0.22 


0.25 


-0.05 


0.00 




63.66 


1.60 


-1.52 


-0.93 


118.17 


0.66 


-0.18 


0.04 


0.22 


0.25 


-0.05 


0.00 




63.23 


1.61 


-1.48 


-0.93 


118.15 


0.66 


-0.18 


0.04 


0.22 


0.26 


-0.05 


0.00 


10 


59.95 


3.65 


-3.04 


-2.06 


102.61 


1.16 


-0.68 


0.16 


0.85 


0.71 


-0.18 


0.01 




60.10 


3.64 


-3.08 


-2.05 


102.62 


1.15 


-0.68 


0.16 


0.85 


0.71 


-0.18 


0.01 




59.99 


3.61 


-3.11 


-2.04 


102.59 


1.13 


-0.67 


0.16 


0.85 


0.71 


-0.18 


0.01 




59.46 


3.64 


-3.04 


-2.05 


102.59 


1.14 


-0.67 


0.16 


0.85 


0.72 


-0.18 


0.01 


25 


50.90 


8.13 


-6.31 


-4.88 


80.63 


1.79 


-2.80 


0.68 


3.71 


2.56 


-0.76 


0.09 




51.04 


8.16 


-6.42 


-4.86 


80.59 


1.77 


-2.80 


0.69 


3.72 


2.57 


-0.75 


0.09 




50.88 


8.09 


-6.51 


-4.84 


80.56 


1.73 


-2.80 


0.68 


3.72 


2.56 


-0.75 


0.09 




50.41 


8.24 


-6.37 


-4.85 


80.63 


1.74 


-2.75 


0.67 


3.73 


2.61 


-0.75 


0.09 


50 


40.54 


10.70 


-9.67 


-8.25 


62.77 


2.11 


-6.43 


1.73 


8.97 


5.89 


-1.63 


0.30 




40.56 


10.81 


-9.80 


-8.25 


62.64 


2.09 


-6.45 


1.72 


8.98 


5.88 


-1.61 


0.31 




40.35 


10.76 


-9.96 


-8.24 


62.62 


2.00 


-6.45 


1.72 


8.97 


5.87 


-1.62 


0.31 




40.18 


11.13 


-9.89 


-8.26 


62.78 


2.03 


-6.31 


1.69 


9.00 


6.00 


-1.60 


0.30 


100 


26.78 


8.46 


-14.52 


-13.24 


43.23 


2.42 


-12.24 


3.90 


17.27 


10.94 


-2.58 


0.76 




26.44 


8.57 


-14.42 


-13.30 


42.98 


2.44 


-12.26 


3.83 


17.26 


10.89 


-2.54 


0.76 




26.18 


8.65 


-14.59 


-13.33 


42.95 


2.25 


-12.31 


3.85 


17.22 


10.88 


-2.54 


0.77 




26.32 


9.22 


-14.91 


-13.38 


43.18 


2.36 


-12.07 


3.79 


17.12 


11.21 


-2.54 


0.74 


150 


16.93 


3.69 


-18.65 


-17.46 


30.72 


2.75 


-16.49 


5.79 


22.12 


13.84 


-2.80 


1.12 




16.33 


3.67 


-18.23 


-17.56 


30.47 


2.83 


-16.45 


5.77 


22.15 


13.80 


-2.79 


1.11 




16.02 


3.83 


-18.32 


-17.59 


30.40 


2.59 


-16.61 


5.78 


22.05 


13.78 


-2.78 


1.12 




16.11 


4.47 


-18.91 


-17.68 


30.67 


2.82 


-16.30 


5.68 


21.72 


14.27 


-2.82 


1.08 


200 


8.94 


-1.44 


-22.18 


-21.30 


21.22 


3.13 


-19.71 


7.29 


24.50 


15.46 


-2.70 


1.33 




8.27 


-1.60 


-21.51 


-21.41 


21.08 


3.27 


-19.62 


7.43 


24.61 


15.48 


-2.73 


1.28 




7.92 


-1.42 


-21.52 


-21.38 


20.98 


3.03 


-19.85 


7.38 


24.49 


15.42 


-2.71 


1.31 




7.83 


-0.74 


-22.15 


-21.55 


21.31 


3.40 


-19.46 


7.26 


23.96 


16.01 


-2.79 


1.27 


250 


1.96 


-6.51 


-25.14 


-24.84 


13.39 


3.56 


-22.21 


8.53 


25.40 


16.39 


-2.49 


1.35 




1.48 


-6.69 


-24.28 


-24.93 


13.46 


3.70 


-22.13 


8.78 


25.48 


16.47 


-2.56 


1.26 




1.10 


-6.55 


-24.22 


-24.78 


13.35 


3.55 


-22.34 


8.68 


25.46 


16.39 


-2.54 


1.32 




0.86 


-5.80 


-24.71 


-25.08 


13.81 


4.05 


-21.84 


8.55 


24.81 


17.02 


-2.61 


1.28 


300 


-4.46 


-11.47 


-27.58 


-28.07 


6.60 


4.03 


-24.14 


9.69 


25.45 


16.95 


-2.30 


1.19 




-4.43 


-11.46 


-26.52 


-28.17 


7.00 


4.10 


-24.21 


9.82 


25.32 


17.08 


-2.36 


1.06 




-4.84 


-11.43 


-26.43 


-27.85 


6.90 


4.12 


-24.26 


9.73 


25.52 


16.99 


-2.35 


1.14 




-5.17 


-10.57 


-26.69 


-28.31 


7.55 


4.74 


-23.64 


9.63 


24.84 


17.64 


-2.36 


1.12 


350 


-10.58 


-16.39 


-29.66 


-30.97 


0.50 


4.57 


-25.57 


10.96 


25.08 


17.31 


-2.18 


0.87 




-9.70 


-15.90 


-28.25 


-31.14 


1.36 


4.47 


-26.00 


10.60 


24.48 


17.44 


-2.19 


0.69 




-10.10 


-16.01 


-28.17 


-30.60 


1.30 


4.74 


-25.74 


10.59 


25.02 


17.40 


-2.18 


0.81 




-10.47 


-15.02 


-28.15 


-31.27 


2.20 


5.45 


-25.00 


10.55 


24.40 


18.03 


-2.09 


0.82 
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TABLE VI. Deuteron properties of the potential models: Z?/S'-ratio rj, asymptotic S'-state 
normalization As in fm^^^^, wave function normalization iV^ in fm~^, D-state probability Pd in 
%, and quadrupole moment Qd in fm^. 

Nijm I Nijm II Reid93 Nijm93 

^ 0.0253 0.0252 0.0251 0.0252 

As 0.8841 0.8845 0.8853 0.8842 

0.7821 0.7828 0.7843 0.7823 

Pd 5.664 5.635 5.699 5.755 

Qd 0.2719 0.2707 0.2703 0.2706 
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FIGURES 



FIG. 1. The mixing parameter ei of the various potentials and of the Nijmegen PWA93. The 
shaded band denotes the statistical error on ei as obtained in the Nijmegen partial-wave analysis. 

FIG. 2. The new Nijmegen potentials in the np singlet partial waves up to J = 3. Solid line: 
Nijm I, dashed line: Nijm II, dotted line: Reid93. 

FIG. 3. The new Nijmegen potentials in the np triplet uncoupled partial waves up to J = 3. 
Solid line: Nijm I, dashed line: Nijm II, dotted line: Reid93. 

FIG. 4. The new Nijmegen potentials in the np triplet coupled partial waves with J = 1 and 
J = 2. Solid line: Nijm I, dashed line: Nijm II, dotted line: Reid93. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9406039vl 



This figure "fig2-l.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9406039vl 



This figure "fig2-2.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9406039vl 



This figure "fig2-3.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9406039vl 



